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£^ ' We investigate non-perturbative pair production from vacuum (the Schwinger 

\f} ' effect) in the focal region of two counter-propagating, ultra-short laser pulses 

.^j- , with sub-cycle structure. We use the quantum kinetic formulation to calculate 

^-j- ■ the momentum spectrum of created particles and show the extreme sensitivity 

to the laser frequency ui, the pulse length t and the carrier-envelope absolute 
phase <f>. We apply this formalism to both fermions and bosons to illustrate the 
influence of quantum statistics in this type of electric background field. 
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' 1. Introduction 

H . 

Non-perturbative electron-positron pair production due to the instability of 
the vacuum in the presence of strong external electric fields — the so-called 
Schwinger effect — has been a long-standing prediction of quantum elec- 
trodynamics (QED) 1-3 but has not been observed yet. This effect was first 
considered for spatially homogeneous and static electric fields. The rate is 
exponentially small, with the scale set by the critical field strength which 
is of the order of E ct = m 2 c 3 /eh « 10 18 V/m. While the production of con- 
stant electric fields of this order is rather unrealistic, recent developments 
in laser technology have raised hopes to approach the Schwinger limit in 
the focal region of colliding laser pulses — either at optical high-intensity 
laser facilities such as ELI or in X-ray free electron laser (XFEL) systems. 
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In this investigation we model the electric field produced in the focal 
region of two counter-propagating laser pulses by assuming that the scale 
of spatial variation of the electric field is much larger than the Compton 
wavelength. Thus, we approximate the experimental situation by a spatially 
homogeneous electric field E(t) = (0, 0, E(t)), represented by an oscillatory 
field with a temporal Gaussian envelope: 



E(t) = E cos(LJt + ^)exp(-^ s S j . (1) 

For this type of electric field there is a simple analytic expression of the time- 
dependent vector potential A(t) — (0,0, A(t)) in terms of complex error 
functions. 5 Due to the appearance of such a variety of physical parameters - 
the field strength Eq, the laser freqency u>, the pulse length parameter r and 
the carrier-envelope absolute phase (carrier phase) - we are faced with a 
rather complicated interplay between various scales, which ultimately leads 
to distinctive signatures in the momentum distribution of produced pairs. 4 ' 5 

2. Quantum Kinetic Equation 

The Schwinger effect is a non-equilibrium, time-dependent quantum pro- 
cess and hence quantum kinetic theory provides an appropriate framework. 
The quantum kinetic formulation arises as a rigorous connection between 
kinetic theory and mean-field approximation to scalar QED (sQED) and 
QED. 6,7 The key quantity in this approach is the momentum distribution 
function f±(k, t) which satisfies a non-Markovian quantum Vlasov equation 
including a source term for particle-antiparticle pair production: 

^/±&*) = HT^ / dt'W ± (t')[l±2f ± (k,t')\co S 2 fdt"w{t") 

— oo 

(2) 

Denoting bosons with (+) and fermions with (— ), W±(t) are given by 

^) = ^f - W-U-H^. ,3, 

with e being the electric charge, k = (k±,k\\) is the canonical three- 
momentum vector and p\\(t) = fcy — eA(t) is the kinetic momentum along 
the electric field direction. e\ = m 2 + k\ is the transverse energy squared 
and to 2 {t) = + P 2 (t) characterizes the total energy squared. It is abso- 
lutely crucial to note that f± (k, t) has physical meaning as the distribution 
function of real particles only at asymptotic times t — > ±oo. 
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3. Quantum Statistics Effect 



We consider the subcritical field strength regime Eq = 0.1E cr , with 
r = 2 • 10 eV -1 , first concentrating on the case of vanishing carrier phase 
<j> = 0. It has been shown in a previous publication 5 that the momentum 
distribution function /_ (k, t) in QED exhibits a distinctive oscillatory struc- 
ture for a = ujt > 4, with the oscillation scale set by the laser frequency 
u. An analogous calculation in the framework of sQED gives a very simi- 
lar result. However, due to the difference in quantum statistics, f-(k, oo) 
shows a local maximum at momentum values at which /+(/c,oo) shows a 
local minimum, and vice versa, as shown in Fig. 1. 

The Schwinger effect in the electric field (1) without carrier-phase <j) has 
been investigated previously in the framework of a WKB approximation, 
together with a Gaussian approximation for the momentum distribution: 8 



d 3 V 
dk 3 



exp 



1 



1 



T kf + k 2 ± 



(4) 



with 7 = muj/eE being the Keldysh parameter. In fact, this approxima- 
tion is too crude in several aspects: First, it does not see the distinctive 
oscillatory structure found in the exact (numerical) treatment; second, the 
Gaussian shape is somewhat broader than the true distribution function. In 
order to explain this discrepancy, we apply the quantum mechanical WKB 
instanton method, 9 for which the momentum distribution is 



'dk 3 



exp 



;-2%) 



with 2S^ = i 



+ k 2 ± + - eA(t)} 2 dt , (5) 



with r being the contour around the branch cut. After a change of vari- 
able, from t to T = —A(t)/Eo we expand the instanton action and obtain 
an approximate solution in terms of an infinite series in powers of the di- 
mensionless variables e = e±/(eEoT) and k = kit/(eEor): 
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1, and the next three terms being given by 
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Fig. 1. Carrier phase <f> = 0; Left: Momentum distribution function f—(k, oo) for kj_ = 
(solid line) in comparison with the WKB-Gaussian approximation Eq. (4) (dashed- 
dotted line) and the WKB instanton result Eq. (5) (dashed line). Right: Momentum 
distribution functions /±(fc,oo) for kj_ = in more detail (solid lines). At momentum 
values, where QED predicts a local maximum, sQED predicts a local minimum and vice 
versa. Note that the WKB instanton result Eq. (5) (dashed line) lies almost perfectly on 
the intersection points, while the WKB-Gaussian approximation (4) does not. 

It is straightforward to calculate even higher order terms. But taking the 
first four terms into account, Eq. (5) already agrees very well with the av- 
eraged envelope of the exact momentum distribution. We point out that 
this averaged envelope lies almost perfectly on the intersection points be- 
tween the momentum distribution of scalar and spinor particles, as shown in 
Fig. 1. Note that none of the semiclassical estimates predicts the oscillatory 
structure in the momentum distribution. 

Finally, considering a carrier phase shift cj) = —ir/2, there are momentum 
values at which no particles are expected to be produced, 5 which is due to a 
resonance phenomenon in the equivalent scattering picture. 10 We emphasize 
that this behavior can be expected for any time-antisymmetric electric field, 
corresponding to a time-symmetric vector potential A(t). Again, comparing 
the results for scalar and spinor particles, we observe that f-(k, oo) shows 
a local maximum at momentum values at which f+(k, oo) shows a local 
minimum, and vice versa, as shown in Fig. 2. 

4. Summary 

The momentum distribution of produced particles is extremely sensitive to 
the physical parameters of a laser pulse. The same qualitative behavior is 
obtained for both scalar and spinor particles, but, due to quantum statistics, 
the oscillatory structure is interchanged. This new effect is complementary 
to the interchange of statistics found in the analogue-thermal distribution 
properties of the QED effective action in electric field backgrounds. 11 ' 12 
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Fig. 2. Carrier phase 4> = — 7r/2; Left: Momentum distribution function f—(k, oo) for 
k± = 0. Right: Comparison of the momentum distribution functions f-(k, oo) (solid 
line) and f+(k, oo) (dashed line), in more detail. 



Here we have considered two identical, colliding laser pulses. Recent 
proposals consider more complicated situations to overcome the strong sup- 
pression of the Schwinger effect. 13-15 Applying the quantum kinetic formu- 
lation, it would be possible to determine not only the total rate but also 
the momentum distribution. 

Acknowledgments: We acknowledge support from the DOC program of 
the Austrian Academy of Sciences and from the FWF doctoral program 
DK-W1203 (FH), from the US DOE grant DE-FG02-92ER40716 (GD), 
and from the DFG grant Gi328/5-l and SFB-TR18 (HG). 

References 

1. F. Sauter, Z. Phys. 69, 742 (1931). 

2. W. Heisenberg and H. Euler, Z. Phys. 98, 714 (1935). 

3. J. Schwinger, Phys. Rev. 82, 664 (1951). 

4. F. Hebenstreit, R. Alkofer and H. Gies, Phys. Rev. D 78, 061701 (2008). 

5. F. Hebenstreit et al., Phys. Rev. Lett 102, 150404 (2009). 

6. Y. Kluger, E. Mottola and J. M. Eisenberg, Phys. Rev. D 58, 125015 (1998). 

7. S. Schmidt et al., Int. J. Mod. Phys. E 7, 709 (1998). 

8. V. S. Popov, JETP Lett. 74, 133 (2001). 

9. S. P. Kim and D. N. Page, Phys. Rev. D 75, 045013 (2007). 

10. C. Dumlu, Phys. Rev. D 79, 065027 (2009). 

11. B. Muller, W. Greiner and J. Rafelski, Phys. Lett. A 63, 181 (1977). 

12. W. Y. Pauchy Hwang and S. P. Kim, Phys. Rev. D 80, 065004 (2009). 

13. R. Schiitzhold, H. Gies and G. V. Dunne, Phys. Rev. Lett. 101, 130404 (2008). 

14. A. Di Piazza et al., arXiv:0906.0726 (2009). 

15. G. V. Dunne, H. Gies and R. Schiitzhold, arXiv.0908.0948 (2009). 



